The effective mass of two-dimensional 3 He 
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We use structural information from diffusion Monte Carlo calculations for two-dimensional 3 He to 
calculate the effective mass. Static effective interactions are constructed from the density- and spin 
structure functions using sumrules. We find that both spin- and density- fluctuations contribute 
about equally to the effective mass. Our results show, in agreement with recent experiments, a 
flattening of the single-particle self-energy with increasing density, which eventually leads to a 
divergent effective mass. 
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PACS numbers: 67.57Pq 

Two-dimensional liquid 3 He is particularly interesting 
because it is, even at zero temperature, not self-bound 
and can, therefore, be studied in a wide density range. 
Although governed by one of the simplest Hamiltonians 
for realistic many-body systems, 3 He exhibits a wide 
range of delicate and complex phenomena which have, 
by-and-large, been resilient to a understanding from the 
underlying Hamiltonian. Only recently, Monte Carlo 
techniques have moved to a point where structuralprop- 
erties have been understood from first principles [ij, |2j . 

Low-energy dynamical properties of 3 He at low tem- 
peratures are phenomenologically described by Landau's 
Fermi-Liquid theory, which establishes relationships be- 
tween observable quantities such as the specific heat, 
the compressibility, and the magnetic susceptibility. Un- 
derstanding the so-called Fermi-Liquid parameters in 
3 He has therefore been a recurring issue in theoretical 
low-temperature research. The calculation of Fermi- 
Liquid parameters in terms of Feynman diagrams is op- 
erationally well defined, but the execution of the theory 
from an underlying microscopic Hamiltonian is far too 
complicated to be practical. Hence, many attempts have 
been made to explain the features of Fermi-Liquid pa- 
rameters within semi-phenomenological models |3J, |5| . 

We examine in this paper physical effects contribut- 
ing to the effective mass in two-dimensional 3 He. This 
work is motivated by a recent sequence of measurements 
that seem to indicate a Mott-Hubbard transition in 
quasi-two-dimensional 3 He atomic monolayers. Techni- 
cally, our calculations correspond to those of Ref. 0, 0, 
but we will use as much information as possible from 
accurate ground state Monte Carlo simulations. 

The relevant quantity for the effective mass is the 
single-particle propagator G(k, uj) in the vicinity of the 
Fermi surface. It is expressed in terms of the proper self- 
energy through the Dyson equation || 



t(k) = h 2 k 2 /2 m is the free single-particle spectrum. The 
physical excitation spectrum is obtained by finding the 
poles of the Green's function in the (fc, w)-plane. Several 
steps are involved in constructing practically useful ex- 
pressions for the proper self-energy E*(fc,w). The first 
step is the derivation of effective interactions. We use for 
that purpose results from diffusion Monte-Carlo calcula- 
tions 0). The structure function can be written as 



S(k) = S p {k) + S <r (k)tr 1 -tr 2 .. 



(2) 



where the components are constructed from the struc- 
ture functions for parallel and antiparallel spins. Both 
quantities are obtained by either directly evaluating the 
expectation value of pkP-k or by Fourier transforming 
the corresponding pair distribution functions. The first 
procedure is more accurate for long wavelength up to the 
size of the simulation box, whereas the latter is appropri- 
ate for medium and short wave lengths. Moreover, one 
can determine the long-wavelength limit of S p (k) from 
the bulk compressibility; we will comment on this below. 
We have taken the Fourier transform of the pair distri- 
bution function for wave numbers k > kp, and direct 
simulation data for long wavelength, and have smoothly 
interpolated these data towards k — > 0. The static struc- 
ture functions at two densities are shown in Fig. ^ The 
density static structure function shows the typical be- 
havior, whereas the spin-structure function depends only 
weakly on the density in the most interesting regime. 

The static structure functions are related to the dy- 
namic response functions through the Too sumrule 
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Above, s £ {p, a} refers to the spin channel. Assuming 
a model like the random phase approximation (RPA), 
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FIG. 1: The density structure functions S(k) (upper curves) 
and correlation part of the spin structure functions Sa-(k) — 
Sp{k) (lower curves) are shown at the densities p — 0.046 A~ 2 
(solid lines) and p = 0.061 A' 2 (dashed line). The zero level 
(dotted line) is included as a guide to the eye. 



we can relate the static structure functions S s (k) 
uniquely to the effective interactions V s (k). The tilde 
in potential indicates that we use dimensionless Fourier 
transforms, V s (k) = p J d 3 rV r s (r)e lkr , correspondingly 
the Lindhard function has the dimension of an inverse 
energy. We note in passing that the random phase ap- 
proximation (0} also satisfies the mi sumrule as an iden- 



tity, whereas the inclusion of at least two-particle-two- 
hole excitations is needed to satisfy higher-order sum- 
rules 

The long-wavelength limit of V p (k) is related to the 
bulk compressibility 



d 9 dE 
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+ V p (0+) . 



(5) 



where E and Ep are the energy per particle of the inter- 
acting and the noninteracting systems, respectively. We 
can thus determine V^(0+) from the equation of state. 

Fig. |21 shows the effective interactions V s (q) defined 
through the relations 10, Q at two representative den- 
sities. These effective interactions resemble the Aldrich- 
Pines pseudopotentials which have been derived in a 
similar spirit. The most prominent features are the same: 
the density-channel interaction is repulsive and can lead 
to a zero sound excitation, whereas the spin-channel in- 
teraction does not. We also note that, similar to the 
spin-structure function, there is relatively little change 
in the spin-channel effective interaction when given in 
dimensionless units. 

Getting back to the self-energy, we assume low-lying 
excitations. In that case, the so-called GOW approxima- 
tion 0, for the self-energy, 



E(fc, E) = u(k) +i T£(2 S + l)J ^M G o (k q E 



should be appropriate. We have split the full self-energy 
into an energy-independent mean field term u(k) and 
the two dynamic, energy dependent portions E' p '(fc,£) 
and E^(k,E) originating from coupling to density- and 
spin-fluctuations, respectively. We have taken for u(k) 
the exchange term of the static density-channel inter- 
action V p (q). One could here in principle also use the 
single-particle spectrum of correlated basis functions the- 
ory, but the basic results are very similar and within the 
limits of the present description. 

The self-energy is conveniently evaluated by Wick ro- 
tation in the complex uj- plane, the salient features have 
been discussed in the literature With the stated 

approximations, one obtains the spectrum 

e{k) = t(k) + E(fc, e(k)) . (7) 

In the numerical applications, we have used the "on-shell 
approximation" e(k) — ► t{k) in the self-energy. It is 
our experience in 3 He- 4 He mixtures that this gives good 
agreement with much more sophisticated implementa- 
tions of the same theory ^jjj. Especially, one might be 



M V?(q) Xs (q, u) = u(k) + (k, E) + SW (k, E) (6) 

I 

led to "dress" the single-particle Green's functions in 
the self-energy © by solving Eq. (J7J. However, single- 
particle Green's functions appear in two locations: One is 
the external propagator spelled out explicitly in Eq. JSJl , 
the other location is the particle-hole propagator. To 
maintain the symmetry between "internal" and "exter- 
nal" propagators, one should apply any modifications to 
either both, or to none. Using a non-trivial spectrum the 
particle-hole propagator is, on the other hand, extremely 
dangerous because one would then violate the mo and 
mi sumrules and, hence, modify the overall importance 
of the dynamic part of the self-energy in an uncontrolled 
way. A systematic improvement of the theoretical frame- 
work is the inclusion of pair excitations |l4| . 

Fig. |3 shows the mean field u{k) and the dynamic 
contributions (k,t(k)) and E^{k,t(k)). All three 
terms are, in the vicinity of the Fermi-momentum, rather 
smooth functions of the single-particle momentum k. 
u(k) has positive slope at the Fermi wave number and, 
hence, decreases the effective mass whereas both dynamic 
contributions cause an effective mass enhancement. The 
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FIG. 2: The density-channel and spin-channel interac- 
tions obtained from the corresponding structure functions 
through the RPA relationship are shown for the densi- 
ties p = 0.046 A" 2 (solid lines) and p = 0.061 A~ 2 (dashed 
lines). The interactions are given in units of h 2 kp/2m. 
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FIG. 4: The figure shows the full on-shell spectrum Q for the 
densities p = 0.026 A" 2 , p = 0.046 A" 2 , and p = 0.061 A" 2 . 
The free single-particle spectrum t(k) is also shown for com- 
parison. All functions have been shifted to be zero at the 
Fermi momentum, all energies are given in units of the Fermi- 
energy of the non-interacting liquid. 



5-shape of the spin-channel term is typical for attractive 
interactions 13 Ell ■ At higher momenta, we found 
similar structures as in three-dimensional 3 He originat- 
ing from a coupling of the single-particle excitation to 
the maxon 0,13; these will not be discussed here. The 
full on-shell self-energy is shown, for several densities, 
in Fig. The most evident feature is the development 
of a saddle around the Fermi wave number, which is due 
to spin-fluctuations and leads, at high densities, to an 
instability. 
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FIG. 3: The figure shows the Fock term tt(fc) (short-dashed 
line), the "density" term E p (fc,i(fc)) and the the "spin" term 
E CT (fe, t(k)) (long-dashed line), and their sum (solid line) of 
the density channel self-energy for the density p — 0.046 A -2 . 
All functions have been shifted to be zero at the Fermi mo- 
mentum, all energies are given in units of Ti 2 k%/2m. 

The effective mass is obtained from the single-particle 
spectrum through 

ft 2 k F _ de(k) 
™* ~ dk k=kp ■ 

Both, experiments |r3. UtI Il8| and our calculations, indi- 



cate that the effective mass increases rapidly with density 
and eventually becomes singular. The primary quan- 
tity that one calculates is the single particle spectrum. 
For a theoretical analysis, it is therefore more conve- 
nient to discuss the inverse, m/m*. Fig. [S] shows, as 
our final result, the density dependence of the effective 
mass ratio m/m* and compares it with the experiments 
of Refs. llTlllSl and|a- The fluctuations of our results are 
due to the statistical uncertainties of the Monte Carlo 
simulations. These are most pronounced in the spin- 
dependent correlations since the spin-structure function, 
S a (k) = Siji(k) — iSjj(/c), is the difference of two quanti- 
ties obtained from DMC simulations, but these fluctua- 
tions do not affect our general result. 

The overall theoretical picture is practically the same 
as the experimental one, although our theory overesti- 
mates the correlation effect somewhat and the instability 
occurs around 0.048 A~ 2 . The experimental ratio m/m* 
is, to reasonable accuracy, a linear function of the den- 
sity which goes through zero between p = 0.051 A -2 and 
p = 0.07 A~ 2 , causing a singular m*. Our theoretical cal- 
culations reproduce the experimental values within about 
10-20 percent when compared with the spectrum of the 
non interacting Fermi gas. This is satisfactory consider- 
ing the simplicity of the G0W approximation |£3l . 

The individual contributions from spin- and density 
fluctuations are also shown in Fig. The mean field 
term and density fluctuations contribute somewhat less 
than spin-fluctuations, but are non-negligible. Both ef- 
fects are individually insufficient to reproduce the ex- 
perimental values. This is the case in both two and 
three dimensions [8j , but the relative importance of spin- 
fluctuations appears to be larger in 2D. 

In both, three and two dimensions, one observes a di- 
vergence of the effective mass at some high density. Ex- 
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though spin-fluctuations are stronger in 2D and we are 
more inclined to associate the effect of density fluctu- 
ations to Feynman-Cohen backflow instead of "localiza- 
tion" . We hesitate to identify the flattening of the single- 
particle spectrum with a Mott transition, it may well in- 
dicate either a transition to "anomalous occupation num- 
bers" or indicate simply a lack of self-consistency of the 
calculation of the spectrum. Quantitative improvement 
will first be sought in a more accurate description of the 
response functions 



FIG. 5: The figure shows the density-dependence of m/m* 
at the Fermi wave number k,F (so lid line). Also shown are 
the results from density-fluctuations (long-dashed line) and 
from spin-fluctuations only (short-dashed line). The dashed 
lines with filled markers show the experimental values of Refs. 
ITU . ITsl and 0. The lines through the theoretical data are 
quadratic fits. 
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trapolating the data of Ref. ElES one observes that the 
divergence of the effective mass would appear in 3D at 
a density 0.03 A -3 , whereas the liquid-solid phase tran- 
sition occurs at 0.023 A~ 3 . The 2D situation is some- 
what different: Ref. implies a divergence of to* at 
about 0.051 A -2 ; earlier data 0, suggest a some- 
what higher density. The film freezes between 0.052 A -2 
@ and 0.063 A" 2 |2l|. 

It seems unlikely that freezing and the singularity of 
the effective mass have the same cause. The divergence 
in to* is due to the increasing importance of spin fluc- 
tuations with density. This is manifested very clearly 
in 2D and also visible in 3D. The theory used here re- 
produces those features, in both 3D and 2D, at a semi- 
quantitative level without the need for phenomenologi- 
cal input. The fact that we obtain a negative slope of 
the single particle spectrum is clearly a consequence of 
the G0W approximation; a self-consistent theory should 
not have solutions for unstable situations. Nevertheless, 
relatively simple approximations have often shown the 
same physics as more sophisticated theories in the stable 
regime, who then simply cease to have solutions beyond 
the point of an instability. We are presently not prepared 
to speculate on "what's beyond" the singularity. 

A second interesting question is why the 2D theory 
apparently overestimates the effective mass, whereas it 
underestimates m* in 3D. One can only speculate that a 
2D model is an obvious simplification of the real physical 
situation of an adsorbed film, and little is known about 
the severity of such an approximation for 3 He. 

We have shown that understanding the value of the 
effective mass in two-dimensional 3 He is - in analogy to 
the more common three-dimensional case - not a simple 
problem, and simple paradigms that try to attribute the 
effect to a single cause are genuinely inadequate. Both 
spin- and density-fluctuations have profound effects, al- 
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